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Abstract 

Brueckner orbitals, and the density of the Brueckner reference-state, are shown to satify the same 
cusp condition - involving the nuclear charges - as natural- and Hartree-Fock-orbitals. Using 
the cusp condition, the density of a determinantal state can be used to determine the external 
potential, if the determinantal state is from either Hartee-Fock or Brueckner-orbital theory, as 
well as, determinant states obtained by many other formalisms that are defined by a one-body 
operator, if a portion of the one-body operator - the portion not associated with the kinetic energy 
or external potential - generates a well behaved function when acting on an occupied orbital. 
Using this relationship involving a determinant and its external potential, a variation of Reference- 
State One-Particle Density-Matrix Theory [arXiv:physics/0308056| is formulated, where the trial 
wavefunctions are universal, in the Kohn-Sham sense, since they do not depend on the external 
potential. The resulting correlation-energy functionals, are also, universal, except for a relatively 
small term involving the portion of the expectation value of the external potential with the trial 
wavefunctions that appears beyond the first order. The same approximate energy functionals that 
were shown to be valid for the previous u-dependent, Reference-State One-Particle Density-Matrix 
Theory [arXiv:physics / 0308084 , are shown to be valid for the current approach, except that the 
use of the LYP and Colle-Salvetti functional appear more natural within the current approach, 
since these functionals are universal ones. And since the BLYP and B3LYP functionals contain 
the LYP functional, these approaches are also better suited with the current approach. 
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I. EXTERNAL POTENTIAL DETERMINED BY THE ONE-PARTICLE DEN- 
SITY MATRIX AND THE PARTICLE DENSITY 



There is a one-to-one correspondence between determinant states and their one-particle 
density-matrices . Because of this correspondence, it is convenient to denote a deter- 
minantal state that is determined by a one-particle density-matrix, say 7, simply by I7). In 
addition, any function, say G, that depends on 7, can be written as G{j); this same notation 
also indicates that G is determined by the corresponding determinant, I7). 

Consider the following noninteracting Hamiltonian: 

n 

^: = E(-i^r. + ^(r.)+^^(x.)), (1) 

i 

where the external potential v is given by a fixed set of point charges 

"(■■) = - E tr^- (2) 

mG{Rnuc} 

and the summation is over the coordinates of the nuclear point charges, denoted by {Rnud; 
furthermore, may be non-local and this operator can depend on the spin-coordinate 
u, where the spatial and spin coordinates are denoted collectively by x; in addition, the 7 
superscript appended to w indicates that this operator may also depend on 7 (or equivalently 

It)). 

Consider a determinantal state, say I7), that satisfies the Schrodinger Eq: 

H]\j) = E,\^), (3) 

where the Hamiltonian Hs is given by Eq. furthermore, the noninteracting eigenstate, 
I7), can be expressed by a unique set of occupied orbitals, denoted by {tpo ^ 7, f-y}', each of 
these orbitals satisfy the following one-particle Schrodinger Eq: 

f^ra^A^) = elr^M, a = a,i3, G ^ 7, /^}, (4) 

where the one-body operator is given by 

f, = -'^V^ + viv) + w'i^). (5) 

In addition, we require the operator to be Hermitian and satisfy 

lim \'R — r\w~'(r,uj)iJjZJr,uj) = 0, for all R. (6) 

IR-rl^O '-rxay V / 



In order to emphasize an exclusive dependence upon R, we modify the hmit in this Eq., 
giving 

hm |r — R|w^(r, ci;)'0^^(r, oj) = 0. for all R. (7) 

r— ♦R 

Let us also mention that the set of unoccupied orbitals - orthogonal to {i/Jo ^ 7, f-y} - is 
denoted by {ipu ^ l^f-y} and, in addition, all of our spin-orbitals ipiai'^) have the following 
form: 

i/jia{yi) = Xia{r)a{u), a = aoY/3, (8) 

where the spatial and spin portions are given by Xio-(r) and cr(a;), respectively, and the spatial 
functions Xjo-(r) are permitted to be unrestricted - two spin orbitals do not, in general, share 
the same spatial function, i.e., {xia 7^ Xip) is permitted. 

Multiplying Eq. (@]) by |r — R| followed by taking a limit of this term vanishing, gives 

^j^jn^nUv;- ^^],.,w^o. (9) 

Y mG{Rnuc} / 

where we have used Eqs. 0, ®5 and (Note that this Eq. is also the cusp condition 
[s,!^, however, in that case, ip].^ is a natural orbital and 7 is the one-particle density-matrix 
of an interacting target-state, say 

In order to obtain an exclusive dependence upon R, we, again, modify the limit in this 
Eq, giving 

lim|r-R| (-|Vr^- Yl l^^-^"^=°- 

\ mGlRnuc} ' 7 

Since the second term vanishes unless (R G {Rnuc}), we have 



lim |r - R| j (-IVr^) - Yl 1 = 0' (^1) 

me{Rnuc} 



which can be written as 



Jim|r-R|^^^(x)-i(-iVr^)^i,(x)= J] S^^^Z^. (12) 

me {Rnuc} 



Defining the left side by 



r(^^,,R) = hm |r - R\r.A^)-' (-iVr^) r.A^), (13) 

r— >R 



we can write 

T{r.^,R)= J2 ^^^-Zm. (14) 

me{Rnuc} 

For a set of spatially restricted orbitals: 

= XlirMu), (15) 

it is readily proven that we have 

T{xl,K)= ^^^^Z^- (16) 

m6{Rnuc} 

Multiplying Eq. by (■(/'^^(x'))*, and summing over all occupied orbitals from the set 
{^o ^ 7, A}' gives 

(-|Vr')- 4h„^J 7(x,x') = 0, (17) 

y me{Rnuc} / 



lim |r — R| 

r^R 



where the one-particle density matrix is given by 



7(x,xo= r.Mir.A^')T (18) 

and it is readily proven that we have 

T(7,R)= Y ^^^A, (19) 

mG{Rnuc} 

where 

T(7,R)= lim|r-R|7(x,x')-^(-|Vr^)7(x,x'). (20) 

r— >R 

Since this expression is invariant to the variable x', we can choose (x' = x), yielding 

T(7, R) = Jim |r - R|7(x, x)-^ ((-| V^^) 7(x, x'))^,^^ • (21) 

Since Eq. (fTT|l is also satisfied by the complex conjugate orbital, '?/'2cr(x)*, it is readily shown 
that we have 

T(7*,R)= Y ^i^H,„^™, (22) 

nie{Rnuc} 



where 

T(7*, R) = Jim |r - R|7(x, x)-i ((-IVr^) 7(x', x))^,^^ . (23) 
Adding together Eq. (fTT^ and (j22I), using (j2II) and (j2Sl), and using the following identity: 
(-iVr^)7(x,x) = ((-|Vr^)7(x,x')),,^,+ ((-iVr^)7(x',x))^,^^, (24) 

we get 

T(p;,R) = 2 J2 ^^^,.Zm, (25) 

me{Rnuc} 

where 

Tip;, R) = Jim |r - R|p;(x)-i {-^V^) p;(x), (26) 
and p^(x) is the spin density, i.e., 

p;(x) = 7(x,x). (27) 

Eq. (fTTjl is also valid for 7(r, a;; r', cu) replacing 7(x, x'); making this substitution and 
summing over the spin- variable u we obtain the same expression, as Eq. (fTTjl . except that 
it involves the spinless density matrix pi, given by 

pi(r,r') = 5^7(r,a;;r',a;), (28) 
and it is readily proven that we have 

r(pi,R) = lr(p,R)= Yl ^HR^^n^, (29) 

mG{Rnuc} 

where T(pi, R) and T(p^, R) are defined by Eqs. (PT|) and ^I^i . respectively; p is the electron 
density, i.e., 

p(r) = pi(r,r). (30) 

Consider the set of (ground and excited) determinantal states, denoted {|7)i,}, that are 
eigenfunctions of Hs, given by Eq. ((T)), where the states from the set, {|7)t,}, are obtained 
from all w~' that satisfy Eq. ((7j), and from all Coulombic external-potentials v, given by 
Eq. 0. From the density of any one of theses states, say p, we can determine its Coulombic 
external-potential v by using Eqs. and ((21). Hence, v is a. unique function of the density. 
In other words, we have f (p), and this function is defined for all densities that are from this 
set of determinantal states, {|7)«}. 
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II. INVARIANCE OF OCCUPIED-ORBITAL TRANSFORMATION 

We now partition the operator into the following four components: 

W = + + <c + (31) 

where the excitation (ex), de-excitation (de), occupied (oc), and unoccupied (un) parts are 
given by the following expressions: 







(32) 




rawer' 


(33) 


"^oc 


= Y '^xa'^La^a', 
waxa' 


(34) 




rasa' 


(35) 



and the occupied- and unoccupied-orbitals are, respectively, given by 

rwa,r.a e {^o^7,/7}> (36) 
rra,rsa ^ H^u ^ 1 , fy} ■ (37) 

The results from the previous Sec. indicate that f is a function of p for any p determined 
from H], given by Eq. - or, equivalently, any p determined from the one-body operator 
f.y given by (0) - when the operator w~' satisfies Eq. ([7j). Using the partitioning method 
given above, Eq. ^ becomes 

lim \r~R\[wl^{r,u) + w2^{r,u)]rw^{r,uj) = 0, for all R, (38) 

r— R|— »0 

indicating that Eq. ((Tj) can be satisfied with any choice of w'^^ and w^^. 
The above relation is satisfied when we have 

\im\r -R\wl^{r,u)'^l^{r,uj) = 0, (39a) 

r^R 

hni |r - R\wl^{r, u;)ijla{^, ^) = 0. (39b) 

It is easily proven that a determinantal state I7) that satisfies Eq. Q - and the corre- 
sponding density p from I7) - does not depend on w"^^; so, when considering the statements 
appearing in the last paragraph of the previous section, we can relax the requirement that 
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Eq. ((Tj) be satisfied, and only require Eq. ()39b|) to be satisfied. In other words, the density of 
a determinantal state, that satisfies Eq. Q, can be used to determine its external potential, 
given by Eq. by using Eq. ^2^ . if ()39b|) is satisfied. An equivalent statement refers to 
the one-body operator f^: The density of a determinantal state can be used to determine 
its external potential, given by Eq. by using Eq. if (l39H|l is satisfied, where the 

orbitals defining the determinantal state I7) are the occupied eigenfunctions of defined 
by Eq. (0). Note that the w^^ and w'^^ portions of the operator w'' are at our disposal, since 
the determinantal state does not depend on these components; ty^g is determined by w2^, 
since is required to be Hermitian. (The Hermitian requirement can be dropped by using 
a biorthogonal basis set.) 



III. HARTREE FOCK DETERMINANTAL STATES 

We now show that the set of Hartree-Fock determinantal states, say {\f}, are members 
of {It)?)}, indicating that their Coulombic external-potentials v can be uniquely determined 
by their electron density, i.e., v{p), by using Eq. (j^ . 

The occupied, canonical Hartree-Fock orbitals satisfy the following single particle Eq: 

AV'L(x) = eL^L(x), a = a,P, (40) 

where the Fock operator is given by 

Ff = -iVr' + ^;(r) + J r,-;f (X2, X2) d^2 + v^i^), (41) 

and the one-particle density-matrix for the Hartree-Fock reference-state has the following 
form: 

f(x,x')= J2 V'L(x)(C.(x'))*; (42) 

furthermore, the exchange operator, v^, is a non-local operator that is defined by its kernel, 
—ri2l- Therefore, for an arbitrary function, say -0, we have 

{);(xi)?/)(xi) = -y dx2rf2^f(xi,X2)V'(x2). (43) 

Equating Eqs. (HO)) and dH)) with (jlj) and (0), for (7 = f), we have 

w^{x) = j r^2^f{x2,X2) dx2 + v^{x), (44) 



and it is easily seen that Eq. (|7j) is satisfied; so, the Hartree-Fock states are members of 
{It)?;}, and we have v{g) where g is the Hartree-Fock density: 

^(r) = X^^(r,^;r,a;). (45) 

The Hartree-Fock Eqs. are usually solved using a iterative, self consistent field (SCF) 
approach, where the (m — l)th iteration is given by 

Ff^-^€7i^) = &ti^), o = a,(5, ^lz^{^o^T.F,^-x} (46) 

and its easily seen that Eq. ((7j) is satisfied for ly""™-!, so all determinantal states determined 
during the SCF approach are also members of {|7)t,}. 

Consider another Hermitian Fock-type operator, say F^, that is given by 

= + vl^ (47) 

where the excitation (ex) portion of the additional exchange-operator i)^^ is zero: 

(^x.)ex = 0- (48) 

Hence, according to Sec. |Hl Eq. ()39b|) remains satisfied and, in addition, the determinantal 
state defined by the occupied orbitals, from F^, is the as same the determinantal state from 
F^; it is a member of {|7)t,}; so, again, the density of this determinantal state can be used 
to determine the external potential Eq. 0, by using Eq. ((221) • Furthermore, since the 
occupied eigenfunctions from F^, given by 

F^^;;(x) = eL^-(x), a = «,/?, ^-e{^o^f,Fa (49) 

differ only by a unitary transformation from the F^ occupied orbitals, {■0o ^ f , Ff }. There- 
fore, and of course, the one-particle density-matrix obtained from these occupied orbitals 
are equivalent: 

f(x,x')= r.M){r..m\ (50) 

where this one-particle density-matrix f is the same one appearing in Eq. ()42|). 
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IV. GENERALIZED FOCK OPERATOR 



Consider a generalized Fock-operator JF^, where its occupied orbitals satisfy 

^i^lA^) = ^L^xa(x), cr = a, /?, V'xa e {^o ^ 7, •5^'} (51) 
and JF^ is given by 

= -|Vr' + i;(r) + j rr2S(x2, X2) rfxs + ^;^(x) + z)-„(x). (52) 
Comparing the above two Eqs. with (@]) and ©, we obtain a new definition for w^: 

w^'{^) = j r^2 1(^2,^2) c?X2 + vl{^),+v2^{:>^). (53) 
and substituting this expression into Eq. ()39b|l . gives 



\\m\r-K\[v2,{r,uJ)l^r^.Ar,i^) = 0, for all R. (54) 

r— »R 

Hence, if this relation is satisfied, the determinantal state defined by the occupied orbitals, 
from JF^, is a member of {|7)v}; so, again, the density of this determinantal state can be 
used to determine the external potential v, Eq. by using Eq. (j^ . 

V. BRUECKNER DETERMINANTAL STATES 

We seek solutions of the time- independent Schrodinger equation, 

H,\'^)=S\^), (55) 

where denotes the Hamiltonian operator defined by the external potential v, where 
the Hamiltonian is independent of the number of electrons when it is expressed in second 
quantization: 

= '^{i(r\h\ja)al^aja + ^ 5Z 5Z {i(^j(^\kcr'l(^')4a^L'(^i'T'(^jT (56) 

iaja iaja ka'la' 

where our Hamiltonian is spin-free; the spin-free integrals are written using chemist's nota- 
tion 0: 

iia\h\ja) = I xl{r)HVl + v{r))x,a{r)dr, (57) 
{iaja\ka'la') = J XiaiT^i)Xja{ri)r\2Xla'i'^2)xia'ir2)dridr2, (58) 
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and the creation and annihilation operators, a]^ and ai^, correspond to the unrestricted 
spin-orbit als, ipia, defined by Eq. (jH)). 

The wavefunction of interest can be generated by a wave operator fl^: 



fi,|7) = (l + X,)|7) = l^), 



(59) 



and the second relation defines the correlation operator, x-y] furthermore, I7) is any deter- 
minantal reference-state that overlaps with the target state: ((7!^) 0). 



Brueckner orbital theory 



IG 
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Q, Q H Q Q, Q, 



20 



2l| is 



a generalization of Hartree-Fock theory that utilizes a single-determinantal state that has 



the maximum overlap with the target state 



23l |. By definition, if |r) is a Brueckner 



reference-state, then the target state, |\E'), contains no singly-excited states from |r): 



(60) 



and the singly-excited states are given by 



(61) 



where the Brueckner-state occupied- and unoccupied-orbitals are, respectively, given by 



e {ipu r}, 



(62) 
(63) 



and this notation indicates that the occupied orbitals determine r; furthermore, the unoc- 
cupied orbitals also determine r since the union of the two orthogonal sets (of orbitals) is a 
complete set. 

Note that, unlike the orbitals that are defined by Eq. (j36p . the occupied orbitals that 
satisfy Eq. ()fi2|l are not completely defined; they are invariant to a unitary transformation; 
similarly, the unoccupied orbitals that satisfy Eq. ()63|) are also invariant to a unitary trans- 
formation. Using a set of these orbitals, the Brueckner one-particle density-matrix is given 
by 



r x,x 



E 

X(tG{V'o^t} 



^L(x) (^L(x'))* 



(64) 
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and, for future use, we mention that the virtual orbitals define the following two-body 
function: 



r(Te{ipu-*T} 



(65) 



where, for a complete set of one-particle functions, the sum of the two gives the Dirac delta 
function: 



5(x, x') = /«r(x, x') + r(x, x') 



(66) 



Since our Hamiltonian, given by Eq. is spin- free, it is easily demonstrated that we 
have 12^ 



(7;;» = 0, for a ^ a' and (7|vl/) ^ ; 



(67) 



hence, we can modify the definition for a Brueckner reference-state, given by Eq. (jHOI), and 
only consider the spin-conserving matrix-elements: 



(68) 



Because of spin symmetry, Eq. ()67|l certainly holds when I7) is a determinantal state that 
is an eigenfunction of the total spin angular- momentum operator, S^, e.g., a closed-shell 
ground-states with spatially restricted spin orbitals. However, this identity should also hold 
in more general cases, since, diagrammatically speaking, the spin state - either a or (3 - 
must be conserved along an oriented path [2J|, and wa and ra' are on the same oriented 
path. In order to simplify our discussions, henceforth, we only consider cases where Eq. fl67|l 
holds; however, the result are easily generalized to the more general case, e.g., when the 
Hamiltonian is spin-dependent. 

Substituting Eqs. and into (jUHjl . sequentially, we obtain 

= (CI*) = {rZ\HM = {rZ\HMr) = (C.|//.|r) + (r^ |if.X.|r), (69) 



and the vanishing of the above matrix elements involving is know as the Brillouin- 



Brueckner condition 
ordered form 24 , 




Writing the operator-product Hj^Xt in normal- 



with respect to the reference state |r), the last matrix element 



of the above Eq. becomes 



{rZ\H.Xr\r) = {rZ\{H,Xr),\r), 



(70) 
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where the the one-body portion, {HyXT)i, can be partitioned in the following manner: 
and where the open (op) portion and remaining (re) portions have the following explicit 



forms 



[(^^X-)i]op = J2^UZala^., (72) 

war a 

rawa rasa waxa 

furthermore, the one-body matrix-elements are defined by 

ui^ = {r,jiH,xr\\r^.), (74) 

and the orbital indices are given by the right side of Eqs. ()62|) and ()63|) : this choice is 
indicated by the r superscripts appended to the summations, i.e., (Note that the 

definition of an open operator given above differs from the definition used by other authors 

■) 
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In the above matrix elements, the ones that do not preserve the spin, i.e., {U^^ for 
a ^ a'), are omitted, since they can easily be shown to vanish for a spin- free Hamiltonian. 
(The vanishing of these matrix elements occurs, diagrammatically speaking, since the spin 



state - either a or (3 - must be conserved along an oriented path 
on the same oriented path.) 

Substituting Eq. ^ into ^ and using ^ and gives 



24| . and ia and ja' are 



(CI (H^Xr), \t) = {rZ\ [iH.Xr),l, |r). (75) 



Since the one-body operator-product [{H^Xt) il^p can also act within the one-body sector of 
the Hilbert space, we have the following identity: 

{rZ\ [iH.Xr),l, \r) = {rj [{H.Xr\l, 1^;.)- (76) 

Substituting Eq. (fTOj) into the Brillouin-Brueckner condition, Eq. (jU^ . and using Eq. 
and (fTBj). and also the following identity: 

{rrJ{Fr)e.\r^J = {rZ\H.m, (77) 
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involving the Fock operator, Eq. ()4ip . yields 

(^;,|(F.)e.|^;j + (^;.| (vi,),^ = o, (78) 

where the introduced correlation potential v^^, by definition, satisfies 

(*cax = [(^.X.)i]op, (79) 

and the operators, (-Fr)ex and {vl^)ex, are defined in an analogous way as w"^^, as indicated 
by Eqs. (jSH) through 

Since the above form of the Brillouin-Brueckner condition, given by Eq. (|78|l . is satisfied 
by all pairs of orbitals involving one unoccupied-orbital and one occupied-orbital, we have 

(j^r) = 0, (80) 

V /ex 

where the generalized, or exact, Fock operator is defined by 

^r = F^ + vl^. (81) 

Comparing this definition of with the one given by Sec. II V| Eq. and using Eq. (jS)), 
we see that, for (7 = r), the two definitions are equivalent, except that in this Sec. we require 
the excitation (ex) portion of the correlation potential v^^ to satisfy Eq. (f?^ : by arbitrarily 
defining the other portions of v^^ we can diagonalize JF^, and this eigenvalue Eq. is given by 
Eq. ^ for (7 = r): 

J^.^;,(x) = £^>;^(x), a = a,/3, r.^e{^o^r,fr}, (82) 

and, furthermore, Eq. ()54|) becomes 

hni |r - R| [vl^{r, uj)l^ ^;,(r, u) = 0, for all R, (83) 

where ('Oco)ex given by Eq. ()79p. Hence, if this relation is satisfied, Bruckner determinantal 
states {|t)} are member of {|7)t,}; so, again, the density of a Brueckner determinantal state 
can be used to determine its external potential v, Eq. by using Eq. (jSHl)- 
Using the results from appendix (^, we have 

hm In - R| [vl^{^^)l^r^.A^l) = lim ki - R| (czK^rM + DZk^i^U^i)) , (84) 

ri~»R ri^R V / 
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where 

V = -^Vr' +i'(ri), (85) 

and there are summations over the repeated indices xa and ra for the orbital sets {ipo ^ 
t,Ft-} and {ipu ^ TjF't}- (The coefficients C^^ and D"^^ are defined by Eqs. ()B38|) and 

jHini).) 

Unfortunately we have been unable to prove that Eq. (jHBI) is an identity by using Eq. (jHlj). 
So, as an alternative approach, consider the case where the above identity, given by Eq. (jEHI), 
is not necessarily satisfied. As in the derivation Eq. (jH)), by multiplying Eq. ()82|) by |ri — R| 
followed by taking a limit of this term vanishing, gives the following identity that must be 
satisfied: 

hm |R- ri|/i,,^;,(xi) + hm In - R| [{;,^„(ri, o;)]^^ ^'^.(ri, ^) = 0, (86) 

ri^R ri— >R 

where we have used Eqs. (1^ . (|H^. and PT|) and, also, omitted the Coulomb and exchange 
terms, since these terms vanish; furthermore, we have used the decomposition of v^^ as 
defined by Eq. (pTTjl . and have chosen ('Oco)oc to be zero, since, according to the discussion 
within Sec. m this portion is at our disposal; the one-particle density-matrix r is invariant 
to this choice. 

Substituting Eq. ^ into dHHI), gives 

hrn^ In - R| (czk,r.M + DZk^rrA^)) = 0, (87) 

where 

C^wl = ^wa,xa + (88) 

Since both terms from the above identity are independent, apparently, we must have 

lim \K-r\K,rxM = 0, (89) 

|r— R|^0 

lim |R-r|/i^iV;^(x) = 0. (90) 

|r— R|^0 

Substituting these relations into Eq. (jH^ proves that Eq. (jH^ is an identity. Hence, the 
density of the Brueckner-determinantal state, |r), can be used to determine the external 
potential v,Eq. (0), by using Eq. Note that Eq. ^ is identical to Eq. ©, the cusp 
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condition, except that the orbitals are now Brueckner; Eq. (|H^ can also be used to prove 
all relations within Sec jl] that appear after Eq. Q, including the one above that states that 
the density of the Brueckner-determinantal state can be used to determine the external 
potential. 



VI. ONE-PARTICLE DENSITY-MATRIX THEORY 



A. Variational Brueckner orbital theory 



Reference-state one-particle density-matrix theory |3 



bital theory 



SQBByyHyyyyyt 
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IS 



341 is based on Brueckner or- 



20 



2l[. Unlike many other 



density functional formalisms based on variants of the Konh-Sham method, the correlation 
operator for this approach is non-local. The approach also uses an energy functional that 
depends on the one-particle density-matrix of a reference determinantal-state, and not the 
exact one from the target state, where the energy functional is partitioned into the exact 
exchange-energy and a correlation-energy functional that is non-universal, since this func- 
tional depends on the external potential v. We now modify this formalism to remove - 
for the most part - the dependence of the correlation energy-functional upon the external 
potential. However, an additional term that describes a portion of the potential energy is 
also included that does not have an analog in Kohn-Sham approaches. However, this term 
can be easily treated once the kinetic energy functional is known or, in many cases, this term 
can be neglected, since it is probably quite small. For convenience, we refer to the previous 
works 30|, |3l|, |3J| as being f -dependent, even though we still retain some t'-dependence in 



the correlation-energy functiona 
In this previous work 



30 



31 



s for the current approach under consideration. 

I, we introduced four t^-dependent trial wavefunctions 
- say I^E'iy'j]), where (77 = I, II, III, and IV) - that are defined with respect to an exter- 
nal potential v and a one-particle density-matrix, where 7 is from a single-determinantal 
reference-state, I7). 

The first trial- wavefunction l^f^^'^) is simply the target state of interest, say |\&ivi,), with 
the single excitations removed: 

\K) = (^-Pn)\^.^), (91) 
where I'ifNv) is a ground-state determined by the external potential v and the number of 
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electrons and, furthermore, the spin-conserved projector for the singly-excited states is 

Pn= E E \^Z){7Z\, (92) 

where we assume that Eq. ^7\i holds, but this requirement can easily be dropped by append- 
ing the states |7^'^') the the right side of Eq. (j^ . Note that the F/i subspace is completely 
determined by I7); P7i ^^^o invariant to a unitary transformation of occupied, or virtual, 
orbitals jl^j ]. 

We now require the state I7) to be a member of {|7)i.}, so Eq. (jHH) becomes 

l^v) = (l-A'\)l^iv.'), 7e{|7).}, 7-iV, (93) 

where the prime superscripts appended to the external potentials, i.e., v', emphasizes that 
the potential defining the states, l^^^^'^/) and |\&]vj;'), may differ from the one determined by 
7 (indicated by 7 w). However, we now restrict these potentials to be equivalent, i.e., 
(f' = v), and generate the target state using a wave operator: 

\^^l) = {l-P^,)\^^.,) = {l-P^,)nM, 7e{|7M, 7-^^, (94) 

and since the external potential is a unique function of 7 (or the density p of 7), we can use 
the function t>(7) to express the wave operator as a unique function of 7: 

n^^ = n^ 7e{|7)4> 7-^^, 7^^, (95) 

and for are target state, we have 

1^^) = 1^^,), 7^v, j^N. (96) 

Therefore, we are assuming that the target state is completely determined by 7. This is a 
reasonable assumption, since the Hamiltonian is completely determined by 7, since 7 gives 
the number of electrons, N, and the external potential v. In addition, however, we must also 
make the assumption that the wave operator is, or - at least, in principle - can be uniquely 
defined so that it generates only one exact eigenstate - the ground state - from all 7 that 
have {-y ^ v). This implies, however, two wavefunctions can differ only by a constant, say 
c, if they are obtained from density-matrices that determine the same external potential: 

1^^) = c|^y), if 7 ^ 7' ^ v', and v = v', (97) 
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where the 7 subscript also indicates the normahzation of the target state: 

(71^7) = 1- 



(98) 



By substituting Eqs. and into we obtain a trial wavefunction that can be 
assumed to be determined by 7: 



and we have 



1*^;') = (1 - P?i) ^.h) = (1 - P^i) 1^.), 7 e {|7).}, 



I^W) = c\w^l), if 7 V, 7' -> v', and v = v', 



(7in') 



(99) 

(100) 
(101) 



The second trial- wavefunction l^&'y"') is defined with respect to the target state expressed 
by an exponential ansatz: (|\1/) = e'^^|7)), where 1^1/^^"') is generated by removing the single- 
excitation amplitudes 5*7 from the cluster-operator S: 



,(5^-57) 



|7), 



where, as in Eq. (jHS)), we have 



The third trial- wavefunction can be generated by its wave-operator: 

^7l7) = 



(102) 



(103) 



(104) 



that can be expressed in an exponential form: (Cl^ = e'^'^|7)), where can be written as a 



sum n-body excitations, with the exclusion of a one-body operator: 



(105) 



The wave operator is a solution to the coupled cluster equations [27|, |28|, |29|, |35|, |36|, |37 



38 



39 



41[ with the single excitation portion removed: 



0, 7 t;, 7 



(106) 



op,cn 



where only the open (op) and connected (cn) portions enter into the relation. This expression 
defines the trial functional using Eq. ()104|) and, again, Eq. is satisfied with 

replacing fl^. 
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The fourth trial wavefunctions l^^i^'^') is not considered here, except to mention that it 
is obtained by solving the configuration-interaction equations p, 0, Q] in an approximate 
way, i.e., by neglecting the single-excitation portion. 

All of the trial states share the property that they contain no single excitations, 

i.e., (P7iI^7') = 0)5 ^iid they generate the target state |\E'^) when their reference state 
satisfies (I7) = |r)), where |t) is the determinantal state constructed from occupied Bruckner 
orbitals. In other words, we have 

1*^"^ = l^r). (107) 

Since the target state is a solution of the Schrodinger equation with an exact energy, 
say £mv, we have 



= ^ J' ,:' = E,[^M+£co[iM = {H,)„ (108) 





\Hy\ 


1^7) 


(^7! 


1^7) 



where the first-order energy Ei is given by the expectation value of the Hamiltonian in- 
volving the reference state, {'y\Hv\^), and the correlation energy Sco is defined above as 
(£^jvt) — El); furthermore, the introduced notation (ifi,)^ indicates the expectation value of 
the Hamiltonian involving the target state {'^■y)- 

Using the trial wavefunctions, we can define variational energy-functionals that depend 
on the one-particle density-matrix: 

E,[^,v] = (H,),, = Ei[j,v] + ^^^'[7,^], (109) 

where we use the notation for the expectation value of an operator, say A, given by 

(■4),, = \:lLv (110) 



and the last relation within Eq. ()109|) defines the correlation-energy functionals E^^^ as 
[Ejj — El); furthermore, the first order energy is given by 

Eih,v] = J [-lV^/9i(r,r')]^,^/r + J viv)p{r) dv + Ej[p] + EMI (m) 

where the Coulomb and exchange energies have their usual forms: 

E^lp] = \j jdvidv2r^^p{vi)p{v^), (112) 

-E^i] = \j j dvidv^Ti^ (p?(ri,r2)p?(r2,ri) + pf(ri,r2)pf(r2,ri)) , (113) 
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and the spin-components of the one particle density matrix are given by 

p?(ri,r2) = Yl XL(ri)(xL(r2))*=7(ri,l,r2,l), (114) 
pf(ri,r2) = Yl X^/3(ri)fe(r2))*=7(ri,-l,r2,-l). (115) 

Substituting Eq. pill) into ()109|) gives the following: 

E,b,v] = j [-iVh(x,x')]^,^^dx + j v{v)p{v)dv + EM + E':^iM, (116) 
where the exchange-correlation energy-functionals are defined by 

E^^[^M = EAl] + E^^bM- (117) 

Returning to our energy functionals, Eq. ()109|) . let the functional derivative of these 
functional yield two-body functions that serve as the kernels of exact Fock operators: 

CJ(xi,X2) = -^J^ (118) 
57(x2,xi) 

= 5{i^2 - xi) + v{v2) + j r2-3S(x3,X3) rfxa^ + z/;-(xi,X2), 

where the kernels of the exchange-correlation operators, //^^"(xi, X2), are obtained from the 
exchange- correlation energy-f unct ionals : 

^r(Xl,X2) = -fj^ = - -r2'7(Xl,X2), (119) 

()7(X2,Xi) 57(X2,Xi) 

where the last relation uses Eqs. pi7j) and the identity: 

-4^^ = -rr2S(Xi,X2) = <(Xi,X2), (120) 

d7(x2,xi) 

and the function i;^(xi,X2) is the kernel of the exchange operator, denoted by {)^. 

Using the variation theorem, and by noting the identi ty g iven by Eq. ()107|) . it becomes 
obvious - as in our previous f-dependent approach Isol . 31 1 - that the minimizing of the 
functionals £^^[7,f], subject to the constraint that the one-particle density-matrix comes 
from a single-determinantal state, yields 

S^, = Er,[T,v], (121) 

£,,[t,v] = E^-^[t,v]. (122) 
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where Ef^^ and Sco the electronic-energy and correlation energy arising from the target state, 
defined by Eqs. pOSj) : furthermore, r is the one-particle density- matrix of the Brueckner 
reference-state |r) that determines v. 

r(x,x')= V^-WCM, r^v, (123) 

and the Brueckner orbitals satisfy the following equivalent conditions: 

{iljXCr'Mn,) =0; ^^e{^o-^r}, e{^„^r}, r^v, (124a) 
(i-r)C|'"r = 0, (124b) 

where these orbitals do not depend oi r] - any trial wavefunction gives the same results. 

A unique set of occupied and unoccupied orbitals is obtained by requiring the occupied 
and unoccupied blocks of C,^'' to be diagonal: 

eV;(x) = CC(x), C e {^o - r}, (125a) 
eV;(x) = C^;(x), e ^ r}. (125b) 

Henceforth, the orbitals sets that satisfy Eqs. ()124|1 and ()125|1 are denoted by {ipl ^ r, Qr^} 
and l'?/'^ ^ r, Cl'''}, indicating that they are determined by r and C,^^\ Since theses orbitals, 
and their energies, can, perhaps, depend on 77, it is more precise to denote then by ipl"^ and 
but we suppress the i] superscripts to keep the notation less cluttered. 
Substituting Eq. ()118|) into Eqs. ()125|) gives generalized, canonical Hartree-Fock Eqs: 



V? + v{v,) + j rr2V(x2,X2) cix2 + />,T(xi)) ^[(xi) = C#(xi) 



lV72 

2 



(126) 



VII. TREATMENT OF THE CORRELATION-ENERGY FUNCTIONALS E, 



The Hamiltonian operator Hy can be written in normal-ordered form j2J, |27|, |28|, |29| with 

respect to the reference state I7): 

= E,[j, v] + {v}, + {-iV^}, + {ri^i}^,, (127) 
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where {v}-y, { — ^V^}^, and {ri2}x2 are terms from the external potential, kinetic energy and 
electron-electron interactions: 

= $^(^^l^li^){«laai^}7' (128) 

icrja 

= (129) 

icrja 

{ru}^ = \^Y1 i'^^j(^\f^^'^^')WicTO-L'(^i<^'^j^}"/ + $^(^^1 (^J + ^x) li^){«Ia«ja}7, 

icrja ka'la' icrja 

(130) 

Furthermore, the Coulomb Vj operator satisfies: 

v'j<P{rr) = I r^^p{v2)4>{v,)dv2, (131) 

and the exchange operator t)^ is given by Eq. (jlHj) . 
Substituting Eqs. (IT77I) into (fTnn|) gives 

bM= [7, v] + [7] + U^l' [7] (132) 

where these terms are identified as the potential, kinetic, and electron-electron-interaction 
contributions to the correlation-energy functionals: 

V^^bM = {{v},U = {{v}),, (133) 

n:'b] = {{-ky%hv = ({-^v^}),,, (134) 

t^rW = {{r^2%h, = {{r^^'Uhv^ (135) 

and where we have also introduced a more condensed notation where the vacuum state is 
understood to agree with the trial wavefunction. 
Similarly, substituting Eqs. ()127|) into (jlOSj) gives 

Sco [7, v] = [7, v] + Tea [7] + Uco [7] (136) 

where 

Ko[7,^] = {{v}),, (137) 

^coW = ({-^V^}), (138) 

Ucob] = {{ru'},^),- (139) 
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VIII. APPROXIMATIONS 



If we know the exact correlation energy Sco for some Brueckner one-particle density- 
matrix, say t' with some external potential, say v' , where r' v', then using Eq. (jl22j) . we 
obtain the following reasonable approximation: 

ETb,v] = g,,[T',v%.=^y=,), r'^v', (140) 



where a similar approximation has been used previously in the f-dependent approach 
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and we are are assuming, as we have previously, that this approximation is most 
appropriate for [r] = III). Similarly, a reasonable approximation for the components of the 
correlation-energy functionals are given by the following prescriptions: 

Vrb,v] = V,o[T',V%^=,y=,), t'^v', (141) 

Trb] = T.o[r'h'=,) (142) 
U^b] = U^ob']ir'=.), (143) 

and the three functionals: V^o, 2^co and Uco, are known for a Brueckner one-particle density- 
matrix t' that determines the external potential v', i.e., r' v'. Of course, we have 
many Brueckner one-particle density-matrices r coming from the same external potential v; 
presumable, we have one r from each A^-electron sector of the Hilbert space, for external 
potentials v with nondegenerate ground states. 

The correlation-energy functionals -E^JJ^ dependence on the external potential v comes 
exclusively form the correlated potential-energy-functional V^™K In many cases it is reason- 
able to assume that the this functionals, and the kinetic energy one T^™\ are small, since 
the potential- and kinetic-energy contributions are treated well in first order. Therefore, the 
following approximation seems reasonable: 

E^''b]^U^ob%'=.), (144) 

and this approximation yields a universal functional in the Kohn-Sham sense - the functional 
does not depend on v. However, even the potential energy contribution to the correlation- 
energy functional, V^o', can be viewed - in a more general sense - as being universal, since if 
this functional is known for an arbitrary external potential, than it is known for other cases, 
since the manner in which this functional depends on the external potential is the same for 
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all systems. On the other hand, if we use a model system to approximate the functionals, 
this will certainly not generate exact functionals for real systems, only approximation ones. 

Using the helium atom as a model system where Uco is presumable known, the previous 
approximation becomes 

^ ^co[rJ(.„.=^), (145) 

where Th,, is the Brueckner one-particle density matrix from the helium atom. Assuming the 
Hartree-Fock one-particle density-matrix, say f^^, is approximately equal to the Brueckner 
one, r„e, we have 

^ t/co[rJ(,„.=^), (146) 



and for a closed-shell systems that use spatially-restricted spin-orbitals, given by Eq. ()A1|1 . 
we have 

[pi] ^ ^co[^j(,,„.=p,), (147) 

where pi is the spinless one-particle density matrix given by Eq. ()A2|) . As demon- 

strated in Appendix ^ we can use the well known approximation for U^^ [Qmc] given by the 
CoUe and Salvetti functional jiJ,!^!, giving 

where this functional, denoted by E^l, uses four empirical parameters that are determined 
using data from the helium atom. In addition, it is reaily verified that this approximate 
correlation-energy functional neglects the potential- and kinetic-energy components, Vco and 
Tco, and these terms vanish when the exact one-particle density matrix, say Fi, from the 
target state I^IZ-y), is equal to the reference state one, 7; these two terms are considered to be 
small, or small enough to neglect, when (Fi ^ 71), as in the approach used when deriving 
the CoUe-Salvetti functional. (Note that Vco vanishes if the density from the reference-state 
is the same as the density from the target state, as is the case for the Kohn-Sham method.) 

Using an identical derivation as in the ^;-dependent approach 3^, it is readily verified 
that the well known density-dependent approximation for the Colle-Salvetti functional j3| , 
given by the LYP functional - at least for closed shell ground states - remains valid for the 
current approach: 

EZ[P]^E'M (149) 
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where the density p dependence is associated with the reference state, and not the target 
state. 

An electron gas defined by an constant external potential is not a member of {|7)t,}, as 
in the case for an electron gas with periodic boundary conditions. However, let us assume 
that we can generalize the functional f (p) so that it yields the appropriate constant value, 
say fg, from the constant-density of an electron gas, say pg\ so, we have {v{pg) = Vg). And 
if we denote the one-particle density-matrix of the Brueckner reference state for an electron 
gas by Tg, we get 

^r[7,^] ^ S,4Tg,Vg]^r,=,,.,=v). (150) 

However, since the correlation energy Sco does not depend on on the constant external 
potential Vg, we cannot make the substitution {vg = v), and so we obtain an approximation 
that yields a universal functional: 

^rw ^ e^'N (151) 

where £^^0""'' is the correlation energy of an electron gas, and this approximation is also 
identical to the approximation used in the f-dependent approach 30|, |3l| , but with a slightly 
different interpretation and derivation. Furthermore, in order to include Tg in the set {|7)i;}, 
we only need to require v{p) to vanishes for any constant density: {v{pg) = 0), where we 
consider two external-potentials that differ by a constant to be equivalent. 

Starting with Eq. (jl51|) . except using a uniform electr on g as, the same correlation energy- 



function used in the local density approximation (LDA) 



3 



46| was shown to be valid with the 
however, in contrast to the Kohn- 



f-dependent approach, for closed shell ground states 
Sham approach, the density dependence of the correlation-energy function is only associated 
with the reference state, and not, in addition, the target state. Furthermore, using an identi- 
cal derivation as in the f-dependent approach, it is easily verified that the correlation-energy 
functional from LDA can also be used, as well, as an approximation for Ej,™^ within the cur- 
rent approach under consideration. Furthermore, since the exchange-energy functional in 
the current approach is identical with the one from the v-dependent method, the exchange- 
energy functionals that are valid in the "y-dependent approach are also valid in the current 
approach, including the Dirac excharige-functional, and the augmentation of this functional 



with the Becke exchange correction 



Hence, as in the f-dependent approach, the LDA 
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and the method known as BLYP are also vahd in the current method. (At least for closed 
shell ground states.) Furthermore, it is readily verified that the B3LYP approach 48, 3| - 
that was demonstrated to be a reasonable approximation within the f-dependent approach 
for closed-shell ground-states js^ - remains valid for current approach under consideration. 

While all functionals that have been shown, so far, to be valid approximations for the 
v-dependent approach, are also valid in the current approach, the use of the LYP and CoUe- 
Salvetti functional appear more natural within the current approach under consideration, 
since these functionals are universal ones that do not have a dependence on the external 
potential. And since the BLYP and B3LYP functionals contain the LYP functional, these 
approaches are also better suited with the current approach. 



APPENDIX A: CONNECTION WITH COLLE SALVETTI FUNCTIONAL 

In order to keep the discussion simple, we only consider closed-shell singlet states that 
are well described by a single determinantal-state, where we use spatially-restricted spin- 
orbit als, given by 

i'jA^) = Xjir)a{uj); a = a,(3. (Al) 

By using these orbitals, it is easily demonstrated that the one-particle density-matrix 7 is 
determined by the spinless one, as indicated by the following relation: 

7(xi,X2) = ^Pi(ri,r2)5^i<^2- (A2) 

Hence, any functional of 7 now becomes a functional of pi; So, if we use the Hartree-Fock 
spin-less one-particle density matrix, say ^1, Eq. (jlOSj) . becomes 

Smv = Ei[gi,v]+ Sco[gi,v]. (A3) 

Substituting Eq. ()L36|1 into this expression for (71 = ^1), and neglected the terms Vco and 
Tco, we have 

S^^^ Ei[qi,v] + Uco[qi], (A4) 

and using Eq. we have 

Smv- J [-iV^^i(r,r')]^,^> + j viT)g{T)dr + Ej[g] + E4g1] + U,o[gi]. (A5) 
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Consider the total electron-electron potential energy, given as the expectation value in- 
volving the target state I^E'pi) and the electron-electron repulsion energy operator {ri2}o'- 

({r,M-. - (A6) 

\^ pi I ^ pi) 

where the operator is given by 

{ri2}o = \ Y1 ^i'i(^j(^\kcr'la')alal^,ai„'aj„, (A7) 

icrja kcr'la' 

and the subscript appended to {rjfg^jo indicates normal-ordering with respect to the true 
vacuum state, | ). Using the Hartree-Fock closed-shell reference-state instead, as the 
vacuum-state, it is readily demonstrated that we have 

({^itlok = {{r^2\2hi + EM + EMi] = UcoIqi] + Ej[g] + EMi], (A8) 

where we have use Eq. ()139|) : in addition, {rjfg^jia is given by Eq. p30|) . where the suppressed 
superscript, 7 - the vacuum state - is set to gi. 

It is well known that the total electron-electron potential energy can also be expressed 
using the (diagonal portion of) the two-particle, spinless density-matrix from the target state 

1m- 








{{ri2}ohi = / nirf{T,,r2)dT,dr2, (A9) 



and we can use the approximate expression, involving a two-body function ip, for the two- 
particle spinless-density-matrix - valid for closed-shell systems - that was derived by CoUe 
and Salvetti 

ri^(ri,r2) = ^2(ri,r2) (l + (/^^(ri, ra) - 2(^(ri, rs)) , (AlO) 
and the two-particle density-matrix, from the Hartree-Fock reference-state, is given by 

hi^i, ra) = i^(ri)^(r2) - ^^i(ri, r2)^i(r2, ri), (All) 

where we also have 

j j r^^H^i,^2)drrdY2 = Ej[g]+EMi]- (A12) 
Substituting Eq. ()A10jl into (|X9|l and using ()A12jl . we have 

{{rl2%)^i = EM + EMi] + EZ[~Qi], (A13) 
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where E^l is the CoUe-Salvetti correlation-energy functional [44. l45l| : 



(A14) 



and this functional, after a series of approximations, is developed into one that does depends 
explicitly on qi |45]. (Note that Q2 is determined by ^1, as indicated by Eq. (jAllI) .) 
Comparing this Eqs. ()A8|) and ()A13|) . gives the desired result: 



(A15) 



and from Eq. ()A4|) . we have 



ij> 



(A16) 



in a gree ment with the Colle-Salvetti electronic energy expresion used in their derivation of 



CO 



44| . where Ei[qi,v\ is the Hartree-Fock energy. 



APPENDIX B: DERIVATION OF Eq. p4|> 

Using the occupied and unoccupied Brueckner orbitals, given by Eqs. (j62j) and (|63|). 
respectively, we can express the one-body operator ('Oco)ex following manner: 

(^co)ex = XI SZ^'^'^I (*™)ex \w(y)al„a^„, (Bl) 

wa rcT 

where, as in {Hxt)i, the vanishing terms involving the matrix elements that do not preserve 
the spin state, , i.e., {ra'\ ('O^o)ex 1^^) fo^' 7^ ^-^^ omitted. The matrix elements in the 
above expression can be computed using the kernel from the operator: 

(^^1 (^^;o)cx k^) = / / V';.(x)<r(x,x')V^;,(x') dxdx'; (B2) 

so, if we have the an expression for the matrix element on the left-hand side, that has the 
general form given by the integral on the right-hand side, we should by able to obtain the 
kernel v^°''(x, x') and, therefore, {v^^)exi^l^ay from the following definition: 

(*;o(x))ex ^;.(X) = I <r(x, xOV'i.lx') dx'. (B3) 

In order to obtain an expression for the term on the left side of Eq. ()86p. below we obtain 
the diagrammatic expansion of the open portion of {Hxt)i, where this operator gives the 
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matrix elements {ra'\vl^\wa) using Eqs. ()79|) and ()B2|) . We then use this matrix element to 
determine 'Oco'0i,cr obtain the identity given by Eq. (jHD). 

It is well known that the correlation operator is given by a linked-diagram expansion, 



where all disconnected pieces are open 
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38L |50| . Since Xt does not contain a one-body 



portion, it is easily demonstrated that the open portion of {Hxt)i is connected - all discon- 
nected pieces from Xt are connected by the Hamiltonian H. 

Using the diagrammatic formalism presented in Appendix [0 and elsewhere js^, consider 
the following example of a diagram that contributes to the open portion of {Hxt)i- 



•\ = ^ri«;r(xi,x2)5f(x2,xiO^/';;;^(xi/)7/';;(xi)a|;^a^^, 

V (B4) 

where the repeated indices, wa and re, are summed over and where the two body function 
is given by 

^(X2, Xi/) = ^£7V(X4, X5)F„/t^(x5, Xg) v(x6, ^4)r:^2l^lyHr{'^2, Xg)^^ v(x3, x'^); (B5) 

furthermore, it is understood that there are no integrations over X2 and Xi' on the right 
side of Eq. flB5j) . since these variables are not repeated indices according to the following 
convention: When determining which dummy indices are repeated indices, indices appearing 
within operators are not counted. So, for example, the indices X2 and x'^ appear only once 
in the above Eq, and not two times, since the dummy indices from the Coulombic operator, 
i.e., X2 (and xe) from rgj^ and x'^ (and X4) from not counted. 

The orbitals ipl^^ and ipl^, presented in Eq. ()B4|) . can be any Brueckner orbitals, as defined 
by Eqs. ()(i2|) and (jHSI)- However, for convenience we choose the canonical orbitals that are 
eigenf unctions of JF^, as defined by Eq. (jH2I)- Using these orbitals the Brueckner one-particle 
density-matrix, r, and the orthogonal function, k,-, are given by 

r(x,x')= (B6) 

xa£{\l)o*-T,FT} 

where these functions are also given by Eqs. and 
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According to Eq. (f7^ . the diagram from Eq. ()B4|1 also contributes to ("Ceo) ex- Comparing 
Eqs. flBl|) and ()B4|) we see that the following term: 

F^iK^(xi, X2)^(x2, xi/)v^;^(xiO^;;(xi), 

contributes to the matrix element (r| {v^o)ex 1"^)' furthermore, and diagrammatically speak- 
ing, removing the incoming and outgoing free- lines from the operator given by Eq. ()B4j) . 
yields 




F^iK^(xi,X2)^(x2,Xi/ 



(B8) 



and this diagram contributes to t>,^^''(xi, xi/), the kernel of (i^co)ex' defined by Eq. ()B2|) : 
furthermore, we use the following diagrammatic representation for this two-body function: 



(B9) 



where, in addition, the non-dummy indices xi and xi' in the above two diagrams, by our 
convention, correspond to the vertices of the omitted outgoing and incoming lines, respec- 
tively. Note that we have labeled these indices in the diagrams above; however, we will often 
omit these labels in similar (kernel) diagrams below. 

As a slight variation of the diagram given within Eq. ()B4|) . consider the following diagram: 



(BIO) 




-FrtKr(xi, X2)5'(x2, Xi/)V';(xi/)^;*(xi)aTa^, 



where the additional factor of ^ comes from the diagonal term, given by Eq. ^C2Q^ . As m 
the diagram within Eq. ()B4|1 . this diagram contributes to ('Oco)cx' corresponding diagram 
that contributes to the kernel i'(^™(xi, xi/) can be expressed in two alternative forms: 




-/Vi/t^(Xi,X2)6f(x2,Xi, 



(Bll) 



where the first diagram on the left side replaces the incoming and outgoing omitted-lines 
with dotted lines; this form gives a visual aid in determining the excitations involved and a 
psudo hole-line for the diagonal term to reside on. 
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The diagrams appearing in Eqs. ()B8|1 and ()B11|1 are examples of one-body kernel-diagrams 
where the omitted outgoing free-line is attached at a Fri vertex and the omitted incoming- 
line is attached at the xi/ vertex of a r~^, operator. (In this particular case we have {j = 4), 
according to Eq. ()B5|1 ). Furthermore, note that the Fock operator is acting upon excited 
orbitals, giving the FT-i/tT-(xi, X2) term. Summing over all diagrams of this type, we have 



F^ifi;^(xi,X2)G'(x2,xi/) 



(B12) 






where the two-body function G(x2,xi/) is obtained from the infinite-order expansion, and, 
for brevity, we have only displayed the first four diagrams of the series. Introducing a 
diagrammatic symbol for G'(x2,xi'), we can represent the above expansion in the following 
manner: 



F^iV(xi,X2)G(x2,XiO 



(B13) 



As a slight variation of the diagrams from the series appearing in Eq. ()B12|) . we also have 
diagrams that have the Fock operator acting upon occupied orbitals, e.g.. 



-F^ir(xi,X2)5'(x2,xi/), 



(B14) 



where (yf(x2,xi/) is given by Eq. ()B5|) . Summing over all diagrams of this type, as in 
Eqs. Im2|l and (fRT8|l . we have 



F„r(xi,X2)/(x2,xi/) 



(B15) 



The diagram sums represented by Eqs. ()B13jl and ()B15j) include all diagram where the 
omitted outgoing free-line is attached at a Fri vertex and the omitted incoming-line is at- 
tached at the xi/ vertex of a r~y operator. Two examples where both incoming and outgoing 
omitted-lines are connected to Fock operators are given by the following two diagrams: 



/Vi V (Xi , X2 ) / (X2 , Xi/ ) /Vi' , 



(B16) 
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F^iT(Xi, X2)m(x2, XiO^^i', 



(BIT) 



where, in these examples, we have 



Z(X2,X1/) = — e^V(x4,X5)F„«;^(x5,X6)Kr(x6,X4)r42V63V(x2,X3)«;^(x3,Xi/), (B18) 

m(x2,xi/) = -^£:;:V(x4,X5)F^5fi;^(x5,X6)fi;^(x6,X4)r43V62Vr(x2,X3)r(x3,Xi/), (B19) 



96 

and again, we can sum over all diagrams of these types: 
F^i«;^(xi,X2)L(x2,xi/)F^i/ = 



F^ir(xi,X2)M(x2,xiOF, 



ri' 




(B20) 



(B21) 



The other two cases of interest involve diagrams where both incoming and outgoing 
omitted lines arc attached to the two-body part of the Hamiltonian, e.g., r^j^, and diagrams 
where the omitted incoming free-line is attached at a /Vi' vertex and the omitted outgoing- 
hne is attached at the Xi vertex of a operator. Examples of these two case are given by 
the following two diagrams: 



(B22) 



(B23) 



where, for these examples, we have 



p(xi,xi/) = j^e^ r(x4,X5)F„«;T-(x5,X6)«;T-(x6,X4)r6i,r4i T(xi,X2)F„r(x2,X3)F„r(x3,xi/). 



n(Xi, Xi/) = j^e^ V(X4, X5)F^5«;^(X5, X6)«;^(X6, X4)r63V4iV(Xi, X2)Frt^(x2, X3)r(x3, Xi 



Summing over all diagrams of these types, gives 

P(xi,Xi/) = .^.^...y^..^. 



N 



(B24) 

') 

(B25) 

(B26) 
(B27) 
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Since the diagrams represented in Eq. (IbT3|) . (pT5|l . (lB2n|l . (|B2T|l . (lB26|) . and (lB27|l 
include all possible diagrams that can contribute to t;^°''(xi, xi'), using these expression and 
Eq. ()B9|) . we have the following diagrammatic and algebraic relations: 



I 



+ " \ + ^ . ^.J (B28) 



^ r 1 V 1 ' 

t';;;"(xi, xi/) = F„K^(xi, X2)G'(x2, xi^ + F^ir(xi, X2)/(x2, xi/) + F^iv(xi, X2)L(x2, xi/)^^' 

+ F„r(xi,X2)M(x2,xr)F„' +P(xi,xiO + Ar(xi,xiOF,i'. (B29) 
Substituting this expression into Eq. ()B3|) . and reordering terms, we have 

(^;o(xi))e.V^;.(xi) = P(xi,XiOV^;.(xiO + iV(xi,XiOF.i'^;.(xiO (B30) 

+ F„/«^(xi, X2)G'(x2, Xi/)V^;^(xi/) + F^ir(xi, X2)/(x2, xi/)^,;;,^(xi/) 

+ F„/«^(xi, X2)L(x2, Xi>)F^^>ipl^^{xr) + /;ir(xi, X2)M(x2, Xi/)F„'V^;;,^(xi/) 

If we substitute the above expression into Eq. (|H^. it is easily demonstrated that the 
first two terms from the above expression vanish, e.g., for the first term, we have 

lim |ri - R|P(xi,xi/)V';^(xiO = 0, for all R, (B31) 

where the variable Xi within P(xi, X2) is the independent variable for functions of the general 
form rJ^T(xi,Xi) or rj^^^/€^(xi, Xj), and the dummy indices in these function: Xj and x^, are 
integrated over; furthermore, and in general, all other variable that appear in diagrams 
that contribute to P(xi,xi/) - e.g., xi, X2, X3, X4, X4, and xe for p(xi,xi/) as presented in 
Eq. ()B24|) - are also dummy integration variables. Hence, since the functions rJ^T{xi Xj) and 
rJ^KT-ixijXi) do not and contain a laplacian term - i.e., — |Vrj'r(xi, Xj) or —^Vl^Kr{'^i,^i) 
- or a singularity - i.e., |ri — Rm|^^ ^ the above identity holds. 
Using a similar analysis, we also obtain the following identity: 

lim In - R|iV(xi,xiOPri'V';.(xiO = 0. (B32) 

ri— +R 

Substituting Eq. ()B30|) into (jHHj) . and using the above two identities, we get 
hm In - R| [vMl^i^U^i) = (B33) 

ri— >R 



lim |ri — R 

ri-»R 



F^/trlxi, X2)G'(X2, Xy)^ljl,^{xi>) + F^iT (xi , X2)/(x2 , Xi>)^IjI^^{xi 
+ F^i V(xi, X2)L(X2, Xi/)F„'^;;,^(xi/) + PriT(xi, X2)M(X2, Xi^)Fn'i'l,„{xy) 
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and this expression can be written as 



lim |ri - R| [{);„(xi)]^^ ^/'^^(xi^ 

ri^R 



(B34) 



lim Iri — R| 



ri-»R 



F„r(Xi, X2)v4^^(x2) +F^iK^(Xi, X2)5^^(X2) 



where 



^^<7(X2) = /(x2,Xi/)?/';^(Xi/) + M(x2,Xi/)F„,'0;^(xiO, 
Byoa{'^2) = G'(x2,XiO'0;^(Xi/) + L(x2,Xi/)F„-'0;^(XiO. 



(B35) 
(B36) 



Substituting into Eq. ()B34|) the Fock operator, Eq. (jH)), it is easily seen that the terms 
involving the Coulomb and exchange operator vanish, so we have 



lini |ri - R| [{)^^„(xi)]^^ ^;^(xi 



(B37) 



lim Iri — R| 



ri^R 



/i„ir(xi,X2)A^^(x2) + /i^,i«;^(xi,X2)S^<^(x2) 



where hm is given by Eq. (j85|l . Using Eqs. (j6l|l and (j65|l . we obtain Eq. (j84j) . where 

CZ = irMY (B38) 

DZ = B^^,) (^;,(X2))* (B39) 

and the terms that do not preserve the spin state, e.g., i^^^', omitted, since these terms 
vanish; furthermore, there are summations over the repeated indices xa and ra for the 
orbital sets {i/Jo ^ t, F^} and {ipu ^ t, Ft}, respectively. 



APPENDIX C: DIAGRAMMATIC FORMALISM FOR THE CORRELATION 
ENERGY £co 

In order to keep the notation less cluttered, for this section we use a combined spin-spatial 
notation for spin-orbital indices; for example, ipwo- is now denoted by ipw- 

A diagrammatic expansion for correlation-energy Sco, or for Xt using Lindgren's formalism 



a 

form 



38, 



5G 



), is easier to obtain when all operators involved are written in normal-ordered 
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5l|. For example, the Hamiltonian, given by Eq. (jHSl); can be written as 



H = E,[T] + {Fr} + {r^,'}r, 



(CI) 
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where 



E,[r] = {r\H\T) 

= 'Y^I'^IMM ~^ — — [wxlxtf]) , (C2) 



{Fr} = E 



{alaj}r, (C3) 



{ruh = lY.^ij\kl]{alalaia,}r, (C4) 

ijkl 

and the integrals are now spin- dependent as indicated by the square brackets [■ ■ ■ ] P| . 
Denoting the one-body portion of H by {Ft}, is appropriate, since this term is the Fock- 
operator, except that the second quantized operators are normal-ordered with respect to 
the |r) vacuum state, instead of the true vacuum | ), as in F^. The two body portion of H 
is denoted by {rfg^}^ emphasizing that this operator is determined by the vacuum 

state |r); furthermore, except for the shifted vacuum, the two-body portion of H is r^2J 
when this operator is expressed in second quantization. 

For a perturbative treatment, we partition the Hamiltonian into a zeroth-order Hamilto- 
nian Hq and a perturbation V: 

H = Ho + V, (C5) 
where we require the reference state |r) to be an eigenfunction oi Hq, a one-body operator: 

Ho\r) = Eo\r), (C6) 
Ho = ^eya,taj, (C7) 

and the zeroth-order Hamiltonian is defined by its matrix elements; we choose them by 
requiring the following relation to be satisfied: 

^ij = Cjri = ^ijj (C8a) 

where 

e;. = 0, (C8b) 

= mfom, (C8c) 
e;. = {rr\f:m, (csd) 
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and the one-body operators, /J and are determined by the reference state |t), but the 
dependence of /j and upon |r) is at our disposal; the orbital subspaces are, again, defined 
by Eqs. ^ and 

Using the above choice, our zeroth-order Hamiltonian becomes 



Ho 



E 



(C9) 



W,X^{lpo^T} 

where the appended r superscript indicates that now depends on the reference state |r). 
A linked diagram expansion for Xt and Sco His known to exist for a zeroth-order Hamil- 
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52 



53 



54 
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56[. A diagonal 



tonian that is a diagonal, one-body, operator [2J, 
form for our one-body operator, H^, is obtained when we choose its orbital sets - {ipo^ t} 
and {ipu^ t} ^ to satisfy the following conditions: 



(ClOa) 
(ClOb) 



where we denote these particular sets of orbitals by {V^o^ t, fo} {i^u^ t, /^}, indicating 
that they are uniquely determined by |r) and their one-particle operator, /J or 
Using these orbitals, can be written as 



Hl = Or + ilr, 



(Cll) 



where these terms - Or and ilr - are the occupied and unoccupied portions of Hq - (Hq) 
and (-f^o)un ^ and are given by the following: 



where our partitioning can be written as 

H = m + Vr. 



(C12a) 
(C12b) 



(C13) 



Using the above notation, our zeroth-order Hamiltonian in normal-ordered form can be 
written as 

H^ = Eo[r] + {dr} + Ur, (CM) 
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where Ur is already normal-ordered; the constant term Eq[t] is the zeroth-order energy of 

H^,\r) = Eo[T]\T), (C15) 

and is given by 

Note that the first-order and the correlation energies, Ei[t] and Sco, do not depend the 
zeroth-order energy i?o ['?"]) 

The perturbation V^, defined by Eqs. ()C13|) . can also be written in normal-ordered form: 

Vr = v; + v{ + v;, (cir) 

where, from Eqs. ()C1|) . and ()C14|) . the individual terms are given by the following expres- 
sions: 



V: = E^[t]-Eo[t], 

Vl = {Fr]-{dr]-Ur, 



(C18a) 
(C18b) 
(C18c) 



The one- and two-body parts of H, {Fr} and {ri2}T, are given by Eqs. ()C3|) and ()C4|) . 
respectively. The Goldstone dia gram matic representation of these operators can be written 



in the following manner 
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5G 
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{Fr} 



(C19a) 



(C19b) 



The one-body part of the perturbation V{ is usually represented by a single diagrammatic 
operator. However, for our purposes, it is convenient to use separate diagrammatic operators 
for the three terms on the right side of Eq. (|C18b|) . where {F^.} is presented by Eq. ()C19a|) . 
Since the other two terms are diagonal, it is appropriate to simply represent them as unfilled 
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arrows: 



-{Or} = V 



-Ur 



A . 



(C20a) 



(C20b) 



In contrast, hole- and particle-lines, by themselves, are represented by filled arrows: Vcind A 



As a slight alternative to the usual approach to eva. 
energy Sco and the correlation operator Xt 
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uate t 



50 



le diagrams of the correlation 



52 



54 
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we associate 



an internal hole-line corresponding to a t/;-occupied orbital with a ■0i„(xi)'0^(x2) factor; 
we associate a particle line corresponding to an r- unoccupied orbital with a '0r(x2)'i/'*(^i) 
factor, where Xi and X2 denote the dummy integration variables that arise from the vertices. 
Using this convention, the sole diagram involving the Fock operator Fr from second-order 
perturbation theory can be evaluated in the following manner: 



^rw) ^ J C?Xi(iX2 ^F^ir^(xi,X2)j • F^2rr(x2,xi), (C21) 




where 



'rw 



(C22) 



and the repeated indices - r and w - are summed over; Ft, denotes the Fock operator Fr - 
given by Eq. PT|) - acting upon (xj); the term {Fn ■ ■ ■)■ indicates that Fn exclusively acts 
within the brackets; furthermore, the wth component of the (one-particle) density-matrix r 
is denoted by 

r^(xi,X2) = ^/^„(xi)?/)* (X2); (C23a) 



the rth orthogonal-component of r is denoted by 

rr(xi,X2) = '0r(xi)?/^*(x2) 
where, for a complete set of orbital states, we have 



(C23b) 



5(xi - X2) = ^r^(xi,X2) + ^r^(xi,X2), 



(C24) 
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which is a shorthand notations for 



5(xi - X2) = 6{ri - r2)5a;ia.2- (C25) 

If we remove the top interaction from the diagram given by Eq. fjC21|) . we see that this 
is a first-order diagram that contributes to the one-body portion of the correlation operator 



Xr B, 



38, 



50|. Since the infinite-order sum of all one-body diagrams for Xr must vanish 
for a Bruckner orbital description, this diagram can be omitted from the expansion for the 
correlation energy Sco However, we will still consider it as a simple example to illustrate our 
approach and notation. 

In order to further compress our notation, we use the convention that all repeated dummy 
indices are integrated over and restrict the Fock operator Fr^ to exclusively act upon the 
first variable of any two-body function, i.e., (^^-.^'(xj, Xj)a(xj, Xj) = a'(xj, Xj)FT-,a(xj, Xj)); 
Eq. ()C21|) can then be written as 

^ = (e™)^^F„T:^(xi,X2)F^2r,(x2,xi), (C26a) 
and the other two diagrams from second-order perturbation theory have the following forms: 



7;{^lwsx) Vi2V34V^(xi,X3)7:,(x3,Xi)r^(x2,X4)r,(x4,X2), (C26b) 



--[£. 



-1^-1^-1 



rwsx) 12 



^34 T^(xi, X3)Tr(x3, X2)t^(x2, X4)rs(x4, Xi), (C26c) 



where 



^rwsx ^rw ~^ ^sx' (C27) 



Let us also mention that when determining which dummy indices are repeated indices, it is 
not necessary to count indices appearing within operators. So, for example, the indices xi 
and X2 appear twice in Eq. ()C26a|) . and not three times, since the dummy indices from the 
Fock operators, i.e., F^-i and F^-j, are not counted. 

The diagonal terms arising from the zeroth-order Hamiltonian, given by —{Or} and —Ut, 
and represented by Eqs. (jC20|) . first appear in third order. For example, the following two 
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diagrams can be obtained by inserting —{Or} and —ilr into the diagram on the left side of 

Eq. jnsnii): 



F^iT^(Xi,X2)/;2T^(x2,Xi), 



(C28a) 



r^rir^(xi,X2)Frtr,.(x2,xi' 



(C28b) 



The hole-hne operator {Or} generates an additional hole line when inserted into a diagram 
and, therefore, a factor of —1 is included when diagram ()C28a|l is evaluated, where this 
factor cancels the —1 factor from — e^. Since this type of cancellation always occurs, as an 
alternative, we associate a factor of for {Or} insertions, and treat {Or} vertices as ones 
that do not generate additional hole lines; ilr is associated with a — factor. Keep in mind, 
also, that these operators generate an additional energy-denominator factor, e.g., e^.^, when 
inserted into a diagram. 

The individual diagrams depend, in part, on each of the Ty^ components, given by 
Eq. ()C23a|) . and the orthogonal components r^., given by Eq. ()C23b|) . In addition, each 
diagram depends on the set of orbital energies {e^}, which are at our disposal. In order to 
make each diagram an explicit functional of the one-particle density-matrix r, given by 

r(xi,X2) = ^r^(xi,X2), (C29) 

to 

and its orthogonal component, Kr, given by 

/t^(xi,X2) = ^rr(xi,X2), (C30) 

r 

where Kr depends, explicitly, on r: 

5(xi - X2) = r(xi,X2) + «;^(xi,X2), (C31) 

we choose all occupied orbitals to be degenerate, with energy el; also, we choose all un- 
occupied orbitals to be degenerate, with energy ej^. With these choices, the zeroth-order 
Hamiltonian, given by Eqs. ()C11|) and ()C12|) . becomes 
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and since this operator is invariant to a unitary transformation of occupied or unoccupied 
orbitals, it no longer depends on /J and - any set of orbitals defining r is appropriate - 
so we can write 

H'o=^:Yl + <E ^l^r. (C33) 

It is easily proven that all perturbative orders, except for the zeroth-order, depend on the 
orbital-energy difference Et, given by 

er = e:- el, (C34) 

and not on the individual orbital-energies, and e^. Therefore, we can choose (e^ = 0), 
and so our only parameter is Et- With this choice we have 

= ErN^, (035) 

where Nt is the number operator for the occupied orbitals. 

If e{V'o— >t} 

and it gives the total number of occupied orbitals when acting on a single determinant. In 
the one-particle Hilbert space, this operator is the projector for the occupied subspace - 
spanned by {^o^ t} - or, the one-particle density-matrix operator: 

Nr^ J2 \i^w){i^w\^r. (C37) 
we{ipo—>-T} 

Using the above two expressions, let us generalize the definition of f: 

f= J2 (C38) 

and write the zeroth-order Hamiltonian in a simplified form, given by 

Hl^ErT. (C39) 

By normal-ordering this expression, we have 

Hl^ErN^ + Er{f], (C40) 
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where A^^ is the number of particles within |r), and from Eq. ()C14|) . we get the following 
identities: 



Eq[t] = SrNr, 
{Or} = er{f}, 

Ur = 0; 



(C41) 
(C42) 
(C43) 



furthermore, our zero- and one-body portion of the perturbation, Eqs. ()C18a|) and ()C18b|) . 
have the following modified forms: 



V; = E,[r]-erK, 
V; = {Fr}-6Ar}- 



(C44a) 
(C44b) 



Eq. ()C43|) indicates that the unoccupied operator, ilr, represented by Eq. ()C20b|) . does 
not appear in the expansion of the correlation-energy S^o', {Or}, represented by Eq. ()C20a|) 
and given by er{f}, is associated with a factor of Sr- Each diagram now becomes an explicit 
functional of r and Kr- For example, the second-order diagrams can be written in the 
following manner: 



^Fr^T{Xi,X2)Fr2Kr{x2, ^1 



(C45) 



-e^ Vi2V34V(xi, X3) v(x3, Xi)r(x2, X4) V(X4, X2), 



(C46) 
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^'^12^^3/'r(^l' X3)/«r(x3, X2)t(x2 , X4) /t^ (X4, Xi), 



(C47) 



where Kr is given by Eq. ((CSH)- Higher order examples are presented elsewhere jSO] and, in 
addition, a method that yields diagrams for the correlation-energy £co that explicitly depend 
on the one particle density-matrix, r. 

It is well known that the correlation operator Yr is given by a linked-diagram expansion. 



where all disconnected pieces are open 
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oratory, 



Using our approach here, these diagrams 
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can be evaluated in an identical manner as the diagrams for the correlation energy 8co] 
Eq. ()C48|) gives an example of a fifth-order one-body Xt diagram: 

= £;:^/Vi«:r(xi,X2)5f(x2,Xi/)?/);(xi/)^;*(xi)aJa^, (C48) 



where the repeated indices, w and r, are summed over and where the two body function is 
given by Eq. ()B5|) . (This diagram is identical to the {Hxt)i diagram appearing in Eq. ()B4|) . 
but that diagram is evaluated slighly different, since there is no energy denominator associ- 
ated with H.) 
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